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ABSTRACT 


The  problem  of  inferring  the  Earth  conductivity  as  a  function  of 
depth  from  surface  measurements  of  the  magnetic  field  is  considered.  Two 
versions  of  a  method  analogous  to  that  of  Bailey  (1970)  are  presented.  The 
main  feature  of  the  method  consists  in  the  fact  that  the  conductivity  is 
given  explicitly  in  terms  to  two  auxiliary  sequences  of  functions  which 
are  found  by  integrating  a  set  of  first  order  nonlinear  differential 
equations. 
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1.  Introduction 

The  inverse  geomagnetic  problem  has  come  to  refer  to  the  problem  of 
retrieving  the  Earth  conductivity  from  temporal  measurements  of  the  magnetic 
field  at  the  Earth  surface.  For  the  case  in  which  the  conductivity  can  be 
assumed  to  be  solely  a  function  of  depth,  Bailey  (1970)  and  Weidelt  (1972) 
have  provided  explicit  solutions  of  the  governing  mathematical  problem.  Their 
results  are  remarkable  in  view  of  the  fact  that  this  inverse  problem  is  non¬ 
linear. 

The  mathematical  formulation  of  the  inverse  geomagnetic  problem  leads 
to  an  inverse  Stt  t-Liouville  problem:  consequently,  any  advance  on  one  front 
can  be  carried  over  to  the  other.  As  a  matter  of  fact,  the  paper  of  Weidelt 
(1972)  is  an  ingenious  adaptation  of  a  method  proposed  by  Gel fond  $  Levitan 
(1956)  for  solving  the  canonical  Sturm-Liouville  problem.  The  present  paper 
falls  in  the  same  category,  and  is  concerned  with  the  adaptation  of  a  method 
for  solving  the  Stuim-Liouville  which  I  have  proposed  (Barcilon  1982) . 

lhe  method  of  solution  which  I  shall  present,  has  some  features  in 
common  with  that  of  Bailey  (1972).  In  particular,  it  provides  an  explicit, 
closed  form  expression  for  the  conductivity.  Also,  its  implementation  requires 
the  integration  of  a  set  of  first  order  nonlinear  equations.  Finally,  it 
requires  some  a-priori  assumptions  on  the  smoothness  of  the  conductivity. 

2.  The  Geomagnetic  Induction  Problem 

We  shall  review  the  induction  problem  very  succinctly  and  refer  the 
readers  to,  for  example,  Le  Mouel's  (1976)  article  for  an  extensive  discussion 
of  the  topic. 

In  a  nutshell  the  situation  is  a1-  follows:  time  varying  magnetic 
fields  of  external  origin,  induce  electric  fields,  and  hence  electric  currents, 
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inside  the  Earth.  These  electric  currents,  which  are  also  time  varying,  in 
turn  generate  magnetic  fields  of  internal  origin.  Surface  measurements  of 
magnetic  fields  can  be  processed  in  such  a  way  as  to  distinguish  between  the 
externally  and  the  internally  generated  parts.  The  ratios  of  these  fields 
contain  information  about  the  conductivity  of  regions  through  which  the  currents 
have  flowed. 

The  following  simplified  version  of  Maxwell’s  equations  describes  the 
situation: 

B.+  Vx  E  ■  0  ,  (1) 

«*t  •  v 

—  Vx  B  »  oE  ,  (2) 

U0  *  • 

V*B  «  0  .  (3) 

In  view  of  the  characteristic  time  scale  of  the  phenomenon,  we  have  made  the 
quasi-static  approximation  in  (2)  and  neglected  the  displacement  current. 

Also  in  (2),  the  electric  current  is  written  via  Ohm's  law,  o  being  the  con¬ 
ductivity  which  is  solely  a  function  of  radius  r.  The  magnetic  permeability 
has  also  been  assumed  to  be  a  constant  uQ,  the  free  space  permeability. 

Because  of  the  solenoidal  nature  of  the  magnetic  field  J3,  it  is  well 
known  that  (see  e.g.  Stern  1976) 

B  ■  -  Vx  Vx  p  r  -  Vx  q  r  (4) 

■»  «•  «• 

where  p  (r,9,<j>,t)  and  q(r,0,<|i,t)  are  the  poloidal  and  toroidal  scalar 
potentials.  Only  the  poloidal  potential  enters  into  the  induction  problem. 
Substituting  the  above  expression  for  8^  in  terms  of  p  only  in  (1)  and  (2) 
and  eliminating  E,  we  get: 


p  ■  m0<j  pt  ,  for  r  <  a  , 


a  being  the  Earth  radius,  whereas  in  the  non-conducting  exterior: 


7*p  ■  0  ,  for  r  >  a  . 


At  the  Earth  surface,  the  three  components  of  B  must  be  continuous  because 
of  our  assumption  about  the  permeability.  As  a  result 


[b^  *  [B01  *  *  0  > 


where  [  ]  stands  for  "jump  across  r«a",  and  the  subscript  r,e,$  refer  to 
the  spherical  coordinates.  Since  o  depends  solely  on  r,  we  write  the 
poloidal  scalar  as  a  series  of  spherical  harmonics  {Y^(0,<fi)}  and  as  a 
Fourier  transform  integral  over  frequency,  namely 


4r  e'iU>t  du  l  I  f?Cr,«)  ,  r  <  a  , 

£■1  m»-£ 


•  ^  »  Xr 

e  dW  l  l  [ 

^  5,-1  m— ft 


r  J  M* 


/a\U1  m 

-V  I)  ]*>.*).  r>a. 


In  the  above  formula,  e  is  the  frequency,  E™(u)  and  I®(o»J  are  for  the  time 
being  coefficients  entering  in  the  most  general  solution  of  Laplace* s  equation 
(6) .  In  order  for  p  to  satisfy  (5)  we  must  require  that 


(9) 


df* 

r2  -jj*  -  4(4+l)f^  ♦  iwour2of^  ■  0  ,  for  r  <  a  . 

By  means  of  this  representation  for  p,  we  can  get  expressions  for 
the  magnetic  field  outside  the  Barth,  vi2. 


+00 


-iwt 


do.  I  I  [lBj(«a)  (±) 
£•1  m—  4  1  V/ 


4- 1 


(£+1)  I 


>>  (?)  'i  A  i 


(10) 


c4+2  3y£  A  3Y™ 

"a?  JL  +  sineH  1 


[E”Cu)(r/  +  I4C“)(f)*  ]  9  *  msr-  *  1  . 


and  inside  the. Earth 


e-iu)t  du  I  (4(4+1)  fj  (r,<o)  Y®  r 


4,m 


3Y® 


3Y, 


F3?rf?(r-“>  hi  i  *  JOT  li  *1 ) 


(ID 


where  r,  0  and  <fr  are  unit  vectors  in  the  r,  9  and  <fr  directions. 

•m  **  «•» 

The  continuity  of  B  across  the  Earth  surface  implies  that 
’(4+1)  f®  (a,«)  -  4E®(u>)  -  (4+1)  lj(u)  , 


(-j^rfHa.u)  *  eJ(u)  +  ij(u) 


(12) 


By  measuring  the  magnetic  field  at  r»a,  we  can  deduce  E®(<o)  and 
I® (w) :  we  shall  look  upon  these  quantities  as  the  data  for  the  inverse 


induction  problem. 
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3.  Properties  of  the  Poloidal  Field. 

The  solution  of  (9)  has  some  properties  which  will  play  a  crucial 
role  in  the  reconstruction  of  o(r).  These  properties  are  associated  with 
the  dependence  of  f^(r,u)  on  the  frequency  w.  In  the  process  of  deriving 
these  properties  we  shall  keep  the  angular  numbers  i  and  m  fixed.  In  fact, 
from  now  on  we  shall  focus  our  attention  on  a  single  pair  (&,m)  and  drop  all 
superscripts  or  superscripts  referring  to  it. 

We  shall  assume  that  the  conductivity  a  is  a  differentiable  function 
of  r.  One  might  be  able  to  relax  this  assumption  at  the  cost  of  complicating 
our  presentation.  For  the  sake  of  simplicity,  we  shall  adopt  this  smoothness 
requirement.  We  should  note  at  this  stage  that  this  assumption  is  weaker 
than  the  analyticity  assumed  by  Bailey.  Incidentally,  in  view  of  this 
assumption  it  follows  that  a(r)  is  bounded,  i.e. 


a(r)  <  M  ,  for  0  <  r  <  a 


It  also  follows  that  the  differential  equation  (9)  has  a  regular  singular 
point  at  r»0  and  hence  it  has  one  solution  which  is  analytic  in  that 


neighborhood  (Coddington  §  Levinson,  1955, 


)  .  It  is  this  analytic 


solution  which  wc  shall  denote  by  f(r,w).  More  specifically,  f(r,w)  is 
that  solution  of  (9)  such  that 


f(r,u)  -  Ar  ,  r  -*•  0 


where  A  is  a  constant. 


We  shall  need  a  better  representation  of  f  in  the  neighborhood  of 
the  origin.  It  is  obtained  by  replacing  a(r)  in  (9)  by  its  value  at 
r=0,  viz.. 
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lff/ 4  t,  u  u 

f(r,w)  -  1.3....(2*+1}A  j  (e  yj  ^  on(0)r]  ,  r  ■*>  0  (IS) 

where  j^(  )  is  the  spherical  Bessel  function  of  order  i  [Abramowitz  4 
Stegun  1972,  p.  437], 

We  next  turn  to  the  dependence  of  f  on  u  and  consider  first  the 
asymptotic  behavior  for  large  w's.  We  resort  to  the  WKBJ  method  (Olver 
1974)  and  look  for  f(r,ui)  in  the  following  form: 


iir/4  k  k  ( 

f(r,u>)  -  exp  [tie  oj 


P  %  ±  al(r) 

dr]  (a*(r)  +  — j- — +  ...) 


Substituting  (16)  in  (9)  and  sorting  terms  according  to  powers  of  w 


we  deduce  that 


+ 


Use  of  the  differentiability  of  a  is  made  at  this  stage. 


a^(r) 


const 

ra^(r) 


(17) 


The  asymptotic  representation  (16)  is  suitable  for  both  solutions  of  (19). 
Therefore,  we  must  select  the  constants  entering  in  (17)  in  such  a  way 
as  to  pick  out  that  solution  which  we  are  calling  f(r,ui) .  To  do  so  we  must 
match  the  asymptotic  representation  for  large  w  with  the  representation  (15) 
for  small  r. 

We  proceed  as  follows:  first  we  note  that  (16)  is  only  valid  for 

r  in  the  interval  (uq  hj  (o)w  ,a) .  This  suggests  that  we  replace  r  in 
-k 

(13)  by  (u0<J(o)uO  o  and  consider  a>>l.  Making  use  of  standard  results 
for  Bessel  functions  (Abramowitz  S  Stegun  1972,  p.  364)  we  deduce  that 
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U  -iTT/4  1  iff/ 4  _ 

f([u0a(oiw]“’o,w)  -  l,3,.,(2W)Ae  a  cost®  a-(fcrl)  y]  .  (18) 


This  implies  that 


t  -u  .w  _w  e±i(A+1)  i 

aj(r)  •  1.3,..(2U1)A  e  *  \  %  *(o)  2 - - - - 


(19) 


Indeed,  the  asymptotic  representation 


f(r,u>)  - 


1.3...(2*+1)A  e 
u^o^(oyro^(r)^~ 


■Lit/ 


4  iff 


cf^fpjdp 


(*♦1)  J  ]  » 


(20) 


agrees  with  (18),  and  hence  (15),  for  small  r's.  We  shall  also  need  the 
asymptotic  form  of  ~  (rf)  which  is  obtained  by  differentiation: 


(rf)  -  -1,3...(2M)A  ff'V)  o^r)  sin[e 


a1*  dr 


(Ul)  i  ]  .  (21) 


We  next  turn  out  attention  to  another  aspect  of  the  same  question  we 
have  been  investigating,  namely  the  dependence  of  f(r,uj)  on  the  frequency. 

We  state  that  f  is  an  entire  function  of  u,  i.e.  an  analytic  function  of  oi 
over  the  entire  complex  ui-plane.  The  proof  follows  a  classical  procedure 
(see  e.g.  Titchmarsh  1962,  p.  6).  We  write  f  as  a  Taylor  series  in  w: 


f(r,u)  *  I  fn<r)  wn 

o 


(22) 


where 


f  "  +  -  f  »  -  f  -  0 

o  r  o  2  o 

r 


8 


and 


also 


fj'.if.'  -ivVr 


n  r  n 


(n>l) 


(23b) 


Si  rml  f (r)  •  A  i  . 
r+o  nv  '  on 


(24) 


Clearly 


and 


f  -  Ar 
o 


n 


iuo 

21+1 


*+2  ^ 

(  Vr'-TT^  °^fn-l(p)  dp 

r  p 


(25) 


It  is  a  simple  matter  to  show  that 


liVr 
o 


l+2n 


I  fn(r)  ■  -  2^4. . . (2n)  •  (21+3)  "(21+5) . . .  (2l+2n+l) 


and  hence  that  the  Taylor  series  converges  everywhere  in  the  complex  en¬ 
plane.  Actually,  we  can  go  further  and  state  that  £(r,u>)  is  an  entire 
function  of  w  of  order  *s+. 


t 


See  (Boas  1954,  p.  8)  for  a  definition  of 


’’order". 
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This  can  be  seen  from  our  previous  discussion  of  the  asymptotic  behavior 
of  f(r,u)).  The  same  can  be  said  about  (rf) ,  i.e.  it  ia  also  an  entire 
function  of  order  *1. 

The  above  results  have  an  important  consequence.  They  enable  us  to 
write  f(r,ij)  and  (rf)  as  infinite  products  (Boas  1954,  p.  22)  of  very 
specific  form,  namely 

cm 

f(r.u)''  •  f(r,o)"Jj~  (1  +  ^(--^y)  ,  (26) 

n»l  n 

at 

(r,ai)  3  “  (rf)  »  iK r,o)  jj  (1  +  Y^-^y)  .  (27) 

n*l  n 

Por  a  fixed  value  of  r,  say  r«o ,  <-iXn(p)}"  and  {  iun (p>  ^  are  the  zeros  of 

f  and  ili  viewed  at  functions  of  w.  In  other  words,  (X  (p)  }?  and  {u  (p)>? 

n  l  n  l 

are  the  eigenvalues  of  the  following  two  eigenvalue  problems: 

u"  +  ~  U'  -  u  ♦  u  X  ou  »  0  ,  re(o,p) 

nrn  *  n  onn 
r 

u  finite  at  r  «  0  ,  (28a) 

u  ■  0  at  r  »  p  , 

and 

V"  +  —  v*  -  v  +  y  v  av  *  0  ,  re(o,p) 

nrn  2  n  onn 

r 

v  finite  at  r  ■  0  ,  (28b) 


(rv)  *  *  0 


at  r  »  p 
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The  reason  ^or  making  this  connection  with  eigenvalue  problems  lies  in  the 
fact  that  {28a)  and  (28b)  enable  us  to  show  that 


An(p)  >  0  ,  for  n»l,2, . . . 


vn(p)  >  0  ,  for  n*l,2, . . . 


Finally,  we  can  easily  find  the  functions  f(r,o)  and  ij»(r,o)  appearing  in  (26) 
and  (27).  Thus,  regardless  of  the  conductivity  profile  o(r),  we  can  say  that 
the  radial  poloidal  field  has  the  following  form: 


n«l  n 


(26') 


00 

where  (xn(r)^  is  an  increasing  sequence  of  positive  numbers.  Similarly 


iHr.w)  »  (1+1)  A  r 


(27') 


4.  An  Explicit  Formula  for  the  Conuuctivity 

There  are  several  ways  of  obtaining  formulas  expressing  a(r)  in 

terms  of  (A  (r) }?  and  (v  (r) }"  We  shall  present  two  formulas, 
n  l  n  l 

Let  us  consider  the  following  contour  integral: 
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m  \  4^  <*“ 

2,1  J  *2Cr.») 


where  r  is  a  circle  in  the  w-plane  of  infinite  rc-dius.  Using  the  calculus 
of  residues,  we  can  deduce  that 


"  2f(r,-iv  )f(r,-iv  ) 

JW  -  l - 7T 1 - ~ 

n»l  ili  (r,-ivn) 


/here  a  dot  stands  for  differentiation  with  respect  to  u.  On  the  other 
hand,  since  the  path  of  integration  is  an  infinitely  large  circle,  we  can 
replace  f  and  <|>  by  their  asymptotic  representations  as  given  in  (20)  and 
(21),  i.e. 


e  iir/4  \  %  J5 

<i  - s - - —  •  cot  [e  /%'aoi5  crdp 

i  •  y  r^afr)^ 

U)  -*«  o  o 


Clearly 


-  (t+l)  j  ]  du 


J(r)  = 


yQr  o(r) 


Equating  these  two  evaluations  of  J(r)  we  conclude  that 


CO 

.2  x  2, 


-2u  rz  l  v‘(r) 
0  n-1  n 


v  (r) 

n  \  x  _ 
TTt  y  k  A. 
y  j  .  3 


y  _ i _ T 


Vr) 

(1' 


vn(r)  ' 

(1-^ 


A  second  formula  for  a(r)  can  be  derived  in  terms  of  (An(r)}^  and  a  third 
auxiliary  spectrum  {<nCr) >"•  To  that  effect,  it  is  convenient  to  introduce 
two  different  combinations  of  f  and 

F(r,ui)  *  r  ^  f(r,u»)  (« 


and 

0(r,to)  =  r*'+*{<|;(r,uO  -  (2.+l)f  (r,u») } 


C 


It  follows  from  these  definitions  that  F  and  9  are  also  entire 
functions  of  u  of  order  H.  Their  product  representations  which  are 


F(r,<u) 


(1 


and 


0(r,aj) 


•  iyoA  S(r)  &> 


Cl  +  ; 


0) 


ikn(r)- 


where 


S(r)  » 


[  p2(-i,+1')a(p)dp  , 


o 


are  obtained  partly  from  the  definitions  (34)  and  (35) ,  partly  from  considering 
the  equations 


and  partly  from  the  eigenvalue  problem: 


W  "  +  -  W  '  -  W  +  u  K  (.p)ou 

n  r  n  2  o  n  n 

T 


w  finite  at  r  »  0  , 


rw'  ■  Jw  »  0  at  r  »  p, 


The  next  step  consists  in  the  discretization  of  (39) .  This  is  accomplished 
by  replacing  a(r)  in  (39)  by  an  expression  involving  sums  of  delta  functions: 


cr(r)  *  l  3.  5(r-r  ) 

j=l  3  3 


In  other  words,  we  consider  model  Earths  made  up  of  N  conducting  spherical 
shells  located  at  r^,  ^  ...  (Parker  1980).  It  is  convenient  to  define 
rQ=0  and  r^^a.  Denoting  by  Fj  (u>)  and  0^  (u>)  the  values  of  F(r,w)  and 
0(r,ai)  at  r»r.j+0  and  substituting  (40)  in  (39)  we  deduce  that 


Fj  +  1  *  Fj  +  24+1  *•  24+1 

Tj 


,24+1  ^  ®j  * 


(41a) 


Vl  ’  ‘  iuoM  C1’  Vl  Fi*l 


(41b) 


Our  goal  once  again  is  to  obtain  a  formula  for  the  conductivity  at  the 
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of  PN+1C“).  namely 


IS 


But,  as  N***  we  can  show  that  r^+a.  Therefore  we  can  approximate  (43) 
as  follows 


and  consequently 
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“o  C  I  ri 


(vv 


X«M 


,(N).aU(N-l)  , 

Wm)\  _jV 


2 

°Ca)  -  hr  * 


szca) .  x  urrr  WlMa) 

0  &4(i+l)  Xl(&)  II  Kk2(a) 


A  similar  formula  can  be  written  for  an  arbitrary  radius: 


■  “o  Ttttir 1 


.«TT 


“k.iMikw 


i  'k2'r> 


5.  Solution  of  the  Inverse  Geomagnetic  Problem. 

So  far  we  have  derived  results  which  hold  for  general  conductivity 

profiles.  We  would  like  to  exploit  these  results  to  retrieve  o(r) . 

Since  the  auxiliary  spectra  {1  }  and  {v  }  on  the  one  hand  or 

n  n 

{Xn>  and  {cc^}  on  the  other,  enter  so  prominently  in  the  formulas  for  o, 
we  express  the  data  (12)  in  terms  of  these  spectra: 


A E(w)-(*+l)I(cu)  =  A£(£+l)  a*TT  U  ♦  jrer  ) 

n-1  nl  1 


E(u))  +  I(oi)  ■  A(£+l)  a 


00 

HT  (1  +  — 2- 

1  I  U  iv  (a 
n»l  n 


■  (£2+£-l)E(u))  +  (£2+2£+2) I (u)  »  -iuQa' C*+1) AS(a)u)TT (1  +  ikM(ay  ) 
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or  better  still 


|AECoO-(JUl)I(u>)|2  -  A2i2(t+l)2a2*  (1 

1 


M 


|B(w)  +  I(«)|2  -  A2(Ul)2a2ilTr(l  -  — § - )  , 

1  Vw 

00  2 

■(&2+fc-l)E(w)+(fc2+2&+2)  I (u)  |  2  -  u  a‘2(U1)A2S2(a)w2T (1  -  -2, - )  . 

V  VC) 


(49) 


In  practice,  approximations  to  the  first  few  values  of  Xn(a),  v  (a)  and 

•cn(u>)  could  be  obtained  by  approximating  the  left  hand  sides  of  (49)  by 
2 

polynomials  in  u  with  real,  positive  zeros.  Note  axso  that  we  can  extract 
the  value  S(a)  of  the  moment  of  o. 

In  order  to  generate  (^(r)},  (v^fr))  and  (<n(r)}  for  arbitrary 
values  of  r,  we  note  that 


f (r,-iAn(r))  *  0  , 

♦(r,-ivn(r))  =  0  , 

9(r,-itcn(r))  *  0  , 

Differentiating  with  respect  to  r  and  making  use  of  the  product  representations 
(26*),  (27*)  and  (36)  we  infer  that 


T, 


dk„  „'r>  k'.'i  Vr)/Vr» 

-a?  *  -V 


(50) 


T  O-  Xn(r)/xk(r)) 

kj*n 
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o» 

f,  ,  „ 

~  <~r~L  ■  Vn(r,ro(r))  JSr- - 


T 

kj*n 


(l-vn(r)/»kCr)) 


I 


and 


d‘n 

IF 


-2(4+1)  o(£) 
S(r) 


KnCr) 


T  (l-<  (r)/X.  (r)) 

k-1  n  k 


"TT  (1"cnCr)/l'kCr)) 

kp<n 


Incidentally,  by  starting  from  the  identity 

F(r,-iXn(r))  »  0  , 

we  could  deduce  similarly  that 


d\>  Vntr) 

dr  2(4+1) 


S(r) 


IT:  l-^_(r)/ic.  (r)) 

k-i  n  k 


Tc  !-*  (r)/X.(r)) 

Wn  n  k 


Two  routes  are  available  to  solve  the  inverse  induction  problem.  The 

first  consists  in  integrating  the  first  order  differential  equations  (50) 

and  (51)  for  {*n(r)}“  and  (vn(r)}"  subject  to  the  initial  conditions 
00  00 

Un(a)>i  ^  deduc®d  from  the  data.  The  conductivity,  which 

incidentally  enters  in  (51),  is  given  by  (33). 

The  second  approach  consists  in  integrating  (52)  and  (53)  and  of 
relying  upon  (47)  for  an  explicit  calculation  of  o(r).  The  fact 


that  S(r)  antars  in  (52)  and  (S3)  is  handled  by  writing 


§-r2(“lV) 

and  recalling  that  S(a)  can  ba  inferred  from  the  data. 

Acknowledgment 

This  research  was  supported  by  the  Office  of  Naval  Research  under 
Contract  N00014-76-C-0034.  Reproduction  in  whole  or  in  part  is  permitted 
for  any  purpose  of  the  United  States  Government. 


(54) 


20 


References 

Abranowitz,  M.  8  Stegun,  I. A.,  1972.  Handbook  of  Mathematical  Functions, 
Dover,  N.Y. 

Bailey,  R.C.,  1970.  Inversion  of  the  geomagnetic  induction  problem,  Proc. 

JR.  Soc.,  135A,  185-194. 

Barcilon,  V.,  1982.  Explicit  solution  of  the  inverse  problem  for  a  vibrating 
string.  J.  Math.  Ansi.  Appl.  (to  appear). 

Boas,  R.P.,  19S4.  Entire  Functions,  Acad.  Press,  N.Y. 

Coddington,  E.A.  8  Levinson,  N.,  195S.  Theory  of  Ordinary  Differential 
Equations,  McGraw-Hill,  N.Y. 

Gelfand,  I.M.  8  Levitan,  B.M.,  1956.  On  the  Determination  of  a  Differential 
Equation  from  its  Spectral  Function,  Amer.  Math.  Soc.  Trans.  1, 

253-304. 

Le-MouBl,  J.-L.,  1976.  L'induction  dans  le  globe,  Chap.  29,  pp.  129-159  in 
Traite'  de  Geophysique  interne,  eds.  Coulumb,  J.  8  Jobert,  G.,  vol.  2, 
Masson  Paris . 

Olver,  F.W.J.,  1374.  Introduction  to  Asymptotics  and  Special  Functions, 
Academic  Press,  N.Y. 

Parker,  R.L.,  1980.  The  Inverse  Problem  of  Electromagnetic  Induction,  JGR, 

85,  4421-4428. 

Stern,  D.P.,  1976.  Representation  of  Magnetic  Fields  in  Space,  Rev.  Geophys. 
Space  Phys.,  14,  199-214. 

Titchmarch,  E.C.,  1969.  Eigenfunction  Expansions,  Vol.  1,  Oxford. 

Weidelt,  P.,  1972,  The  Inverse  Problem  of  Geomagnetic  Induction,  Z.  Geophys., 
38,  257-289. 


